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We have used density functional theory to study the structural stability of surface alloys. Our 
systems consist of a single pseudomorphic layer of M^Ni-^ on the Ru(OOOl) surface, where M = 
Fe or Co, and A*' = Pt, Au, Ag, Cd, or Pb. Several of the combinations studied by us display a 
preference for atomically mixed configurations over phase-segregated forms. We have also performed 
further ab initio calculations to obtain the parameters describing the elastic interactions between 
atoms in the alloy layer, including the effective atomic sizes at the surface. We find that while elastic 
interactions favor alloying for all the systems considered by us, in some cases chemical interactions 
disfavor atomic mixing. We show that a simple criterion (analogous to the Hume-Rothery first law 
for bulk alloys) need not necessarily work for strain-stabilized surface alloys, because of the presence 
of additional elastic contributions to the alloy heat of formation, that will tend to oppose phase 
segregation. 

PACS numbers: 68.35.-p, 68.35.Dv, 68.55.-a 



I. INTRODUCTION 

It has been known since ancient times that alloying two metals can give rise to a new material with properties 
that are improved over those of the constituent metals. For example, alloys can have superior mechanical or 
magnetic properties, an increased resistance to corrosion, or constitute good catalysts. However, not all pairs 
of metals form stable alloy phases. The rules governing alloy formation in the bulk were first formulated by 
Hume-Rothcry.^ The first of these empirical laws states that if the atomic size mismatch is greater than 15%, 
phase segregation is favored over the formation of solid solutions. Thus, many pairs of metals are immiscible 
(or nearly so) in the bulk. 

In recent years, it has become apparent that surface science can extend the chemical phase space available for 
the search for new alloy systems. It has long been known that bulk alloys exhibit surface segregation, so that 
the chemical composition at the surface can differ considerably from that in the bulk. However, the field of 
surface alloying gained additional interest when it was discovered that even metals that are immiscible in the 
bulk can form stable surface alloys as a result of the altered atomic environment at the surfacei^i^ These alloys 
display atomic mixing that is confined to the surface layer or, in some cases, the top few layers. These results 
were explained by Tersoff,^ who argued that in cases where there is a large size mismatch, as a result of which 
the energetics are dominated by strain effects, alloying will be disfavored in the bulk but favored at the surface. 

Subsequently, another class of surface alloys has emerged, where two metals that differ in size are co-deposited 
on a third metal of intermediate size. In such systems, any single-component pseudomorphic layer will be under 
tensile or compressive stress (that may or may not be relieved by the formation of dislocations)'^ however if 
the two elements were to mix, the stress would presumably be considerably relieved. Thus, the strain imposed 
by the presence of the substrate promotes alloying in the surface layer. Some examples of such strain-stabilized 
surface alloys are an Ag-Cu monolayer on Ru(0001)^i^ia Pd-Au/Ru(0001)^ and Pb-Sn/Rh(lll) 

Hitherto, the guiding principle in the search for such systems has been the rule-of-thumb that the (bulk) 
nearest-neighbor (NN) distance of the substrate should be the average of the NN distances of the two overlayer 
elements. However, this simple criterion does not necessarily work. For example, Thayer et al. have studied the 
Co-Ag/Ru(0001) systemJ^iH At first sight, this system would seem to be a good candidate for the formation 
of a strain-stabilized surface alloy, since the NN distance for Ag is larger than that of Ru by 8%, while that of 
Co is smaller by 7%. However, instead of forming an atomically mixed structure, it was found that the stable 
structure consisted of Ag droplets surrounded by Co. After doing a combined experimental and theoretical 
study, these authors concluded that chemical bonding between Ag and Co is disfavored in this system, and the 
observed structure results from a lowering of stress at the boundary between Co and Ag islands. 

In this paper, we examine ten different bi-metallic systems on a Ru(OOOl) substrate. Some of the questions that 
we hope to address include: (i) is it only the mean size of the overlayer atoms that matters, or do individual 
sizes also matter? (ii) can one develop a criterion based on atomic size that will predict whether or not a surface 
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TABLE L Values in A of (qm + a]v)/2, the average of the nearest neighbor (NN) spacings of M and in their bulk 
structures. By comparing these numbers with ag, the NN distance of Ru in the bulk = 2.70 A, one expects Co-Au and 
Co-Ag to form stable alloys, and Pb alloys to be unstable. In this table, all the values used are experimental values, 
taken from Ref. [2ll . 

alloy will form? (iii) how different are atomic sizes at the surface compared to those in the bulk?, and (iv) what 
is the relative importance of elastic and chemical interactions? 

The bi-metallic systems we have considered all consist of one magnetic metal M (Fe or Co) and one non-magnetic 
metal N (Pt, Ag, Au, Cd or Pb), co-deposited on Ru(OOOl) to form a surface alloy of the form M^Ni^x/S, where 
S denotes the Ru substrate. Such systems, involving one magnetic and one non-magnetic element, are of interest 
because alloying can, in some cases, improve magnetic properties. Conversely, in some applications, ^^li^ one 
would prefer that instead of mixing at the atomic level, the system should spontaneously organize into a pattern 
consisting of alternating domains of the magnetic and the non-magnetic element. Ru(OOOl) was chosen as the 
substrate, in part because of its intermediate NN distance, and in part because its hardness and immiscibility 
with the other elements make it less likely that the alloy elements will penetrate into the bulk. The bulk NN 
distances, oa/, of the two magnetic metals Fe and Co, are about 7-8% less than as, the NN separation in the 
Ru substrate, while all five non-magnetic metals we have considered have bulk NN distances, un, larger than 
that of Ru. However, the TV metals we have chosen display a large variation in size: the NN distance in Pt 
is approximately 3% more than that in Ru, while in Pb the discrepancy is 26%. Accordingly, only Fe-Pt and 
Co-Pt fall within the 15% range of the Hume-Rothery criterion for bulk alloys; alloys of Fe and Co with Au 
and Ag fall slightly outside this range, while those with Cd and Pb fall well outside the range. If there is a 
size-dependent trend that determines whether or not alloying is favored, then one might hope that it will show 
up upon examining these ten systems. In Table I, we have given the average NN separation, (a^/ + aN)/2, 
using experimental values for the bulk metals. Upon examining how close these values lie to as = 2.70 A, 
one might expect (using the simple criterion mentioned above) that Fe-Au, Fe-Ag, Co-Au and Co-Ag might 
be good candidates for forming strain-stabilized surface alloys, and Fe-Cd, Co-Cd, Fe-Pt and Co-Pt may be 
possibilities, but Fe-Pb and Co-Pb surface alloys should be highly unlikely to form. As we will show below, 
these simple-minded expectations are not necessarily borne out. 

Of the ten systems we consider in this paper, we are aware of previous studies on only two of them: Co- 
Ag/Ru(0001)fi2iii^ and Fe-Ag/Ru(0001).^^ In both these cases, it was found that chemical interactions dominate 
over elastic ones, and the atomically mixed phase is disfavored. 



II. COMPUTATIONAL DETAILS 

All the calculations are done using ab initio spin polarized density functional theory with the PWscf package of 
the Quantum-ESPRESSO distribution^^ A plane-wave basis set is used with a kinetic energy cutoff of 20 Ry. 
The charge-density cutoff value is taken to be 160 Ry. Ultrasoft pseudopotentials^^ are used to describe the 
interaction between ions and valence electrons. For the exchange correlation functional, a Generalized Gradient 
Approximation (GGA) of the Perdew-Burke-Ernzerhof fornii^ is used. As all the systems are metallic, the 
Methfessel-Paxton smearing technique^^. is used with the smearing width equal to 0.05 Ry. 

Convergence with respect to the basis size and the k-point grid has been carefully verified. For the bulk 
structure calculations, we have used the common crystal phase of each element. The k-points used for Brillouin 
zone integrations form an 8x8x8 Monkhorst-Pack grid^° for bulk calculations, and a 4x4x1 grid for surface 
calculations. To study the surface properties, the supercell approach is used, with a unit cell that includes a slab 
and some vacuum layers. The slab used corresponds to a 2 x 2 surface unit cell, and contains six Ru layers to 
model the substrate. Our results for the energetics were obtained with one alloy overlayer (deposited on one side 
of the substrate) and seven vacuum layers (approximately 17.4 A); we have allowed the alloy overlayer and the 
three topmost layers of Ru to relax, using Hellmann-Feynman forces. However, when performing calculations 



FIG. 1: Top views of the systems studied, for x = (a) 1.00, (b) 0.75, (c) 0.50, (d) 0.25, and (e) 0.00. The green lines 
indicate the boundaries of the 2x2 surface unit cell. S denotes the substrate atoms (gray) and M and A*" denote magnetic 
(blue) and non-magnetic (orange) elements respectively. 



to see how the surface stress of monolayers of M or N on S varied with in-plane distance, the monolayer was 
deposited symmetrically on both sides of the slab, and the central layers of the slab were held fixed, while the 
outer layers on both sides were allowed to relax. 

The (0001) surface of Ru is a closed-packed surface, on which typically one of the hollow sites, either hexagonal- 
closed packed (hep) or face-centered cubic (fee), is energetically preferred. We have allowed for both possibilities. 
The use of a (2x2) unit cell enables us to study five different compositions as shown in Fig. [TJ Because of the 
small size of the unit cell, there is only one distinct configuration corresponding to each composition. 

In this particular study, we have considered only a single, pseudomorphic and ordered layer of an alloy on the 
substrate slab. The general observation that reconstruction in the overlayer occurs only after a certain critical 
thickness of deposited material, validates our assumption that the monolayer (of either a single metal or an 
alloy) remains pseudomorphic. For most systems studied previously consisting of a single metal on Ru(OOOl), 
it is found that the first overlayer of the metal does not reconstruct. However, Ag on Ru(OOOl) is an exception, 
in which a misfit dislocation structure has been observed, even for submonolayer filmsi^ 



III. RESULTS AND DISCUSSION 



The substrate element, Ru, has the hep structure in the bulk. Upon optimizing the geometry for bulk Ru, using 
the experimental c/a ratio of 1.584, we obtain a (which is also the NN distance as) as 2.74 A, which is close to 
the experimental value of 2.70 A.^^ For bulk Fe, Co, Pt, Au, Ag, Cd and Pb we obtain NN distances a = a^^j^. 
of 2.47, 2.49, 2.83, 2.93, 2.95, 3.04 and 3.56 respectively. Again, all these numbers match very well with the 
corresponding experimental values. 

For a single-component monolayer on Ru(OOOl), we find that both the magnetic elements prefer to occupy the 
hep sites; occupying instead the fee sites costs about 75 meV per surface atom. However, for all the non- 
magnetic elements, with the exception of Pt, we find that the fee site is very slightly favored over the hep one, 
with an energy difference of the order of 4 meV per surface atom. In the rest of this paper, we work with the 
structures corresponding to the favored site occupancies for each system. 

Upon depositing the single-component monolayers of either M or N on Ru(OOOl), and relaxing the geometry, 
we find that di2, the interplanar distance between the overlayer and the topmost Ru layer, varies significantly 
depending upon the element constituting the overlayer. Our results for di2 are given in Table HH they may 
be compared with 2.17 A, which is the value of c?^";,, the interplanar distance in bulk Ru. We see that for 
the magnetic elements, di2 < dj^^ik^ whereas for the non-magnetic elements, di2 > dj^^^f^. We also see a similar 
pattern upon examining our results for the surface stress of these systems (see Fig. [2]): the M/Ru systems are 
under tensile stress, whereas all the N/Ru systems are under compressive stress. All these findings are consistent 
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TABLE II: Results for the value of di2, the interplanar distance between the overlayer and topmost Ru layer, for single- 
component monolayers of M or A'^ on the Ru(OOOl) substrate. It is interesting to compare these results with d^^i^. = 
2.17 A, which is the value of the interlayer distance in bulk Ru. 
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FIG. 2: Results for a 



the diagonal component of the surface stress tensor, for a single-component monolayer of 



either M or TV on Ru(OOOl), as a function of the lattice mismatch, defined as (a — as) / { 2^ ). Note for the magnetic 
elements, a < as, and the surface stress is found to be tensile (negative), while for the non-magnetic elements, a > as, 
and there is compressive (positive) surface stress. 



with the idea that the M atoms at the Ru surface would hke to increase their ambient electron density, whereas 
the opposite is true for the N atoms; this is what one would expect from simple size considerations using the 



values of < 



'bulk 



for all the metals. 



For the surface alloys, we find that in every case considered by us, the hep site is favored over the fee site. The 
difference in energy between the two sites varies from 10 to 70 meV per surface atom. Upon relaxing the alloy 
structures, we find that the surface layer can exhibit significant buckling; this follows the trends expected from 
the atomic-size mismatch between the constituent elements. Thus, Pt alloys do not show any visible buckling, 
while Pb alloys show the maximum amount of buckling among all the iV's studied (see Fig. [S]). 

The stability of an alloy phase relative to the phase-segregated phase can be determined by calculating the 
formation energy, AH which is defined as follows: 



AH = Esiab{M,Ni^,/S) - xEsiab{M/S) 
-{l~x)E,iabiN/S), 



(1) 



where Esiab{A) is the ground state energy per surface atom for a single layer of A on the substrate S. When 
AH is negative, the two metals prefer to mix rather than to segregate, and hence the alloy phase is more stable. 

Our results for AH as a function of a composition are presented in Fig. 21 Note that in all cases, we find that 
AH is roughly symmetric about x = 0.5, suggesting that pairwise interactions are dominant. For both Fe and 




FIG. 3: (color online) Relaxed geometries for surface alloys with x = 0.25 for (a) Fe-Pt, (b) Fe-Ag, and (c) Fe-Pb. Here, 
gray, red, purple, green and blue spheres represent Ru, Fe, Pt, Ag and Pb atoms respectively. Note that the amount 
of buckling increases progressively, in keeping with what one would expect upon considering the mismatch between the 
atomic sizes of the constituent elements. 



Co , alloys with Ag are found to be the least stable and alloys with Pb appear to be the most stable. However, 
though Fe and Co have almost the same ag^'^, the values of AiJ, and even the order of stability, are not identical 
in the two cases. Similarly, despite having very close values of a^uik^ ^^'^ display very different behavior: 
alloys of the former are stable, whereas Fe-Ag alloys are right at the boundary of stability, and Co-Ag alloys 
are unstable. These observations support the view that chemical effects may, in some cases, be quite important 
- and even dominate over elastic interactions. Our finding that atomic-level mixing is disfavored for Fe-Ag and 
Co-Ag is in keeping with the observations of previous authors.— li^ii^ We point out that our results underline 
the fact that (cm + cln)/2 « as need not necessarily be a good criterion for atomic-level mixing to be favored 
(see Table HJ. 

It is generally accepted that there are two main contributions to the stability of such surface alloys: an elastic 
contribution, and a chemical contribution.—"^ We would like to separate out the two, if possible. In order to 
do so, we assume that the elastic interactions are given by a sum of NN contributions, with each pairwise term 
taking the form of a Morse potential: 

V,,ir) = AI={1 - exp[-Al^"(r - b^^)]}' (2) 

where r is the distance between the NN atoms i and j, 6'^ is the equilibrium bond length, and Aq and are 
parameters related to the depth and width, respectively, of the potential well. 

For each composition, the elastic energy is written as the sum of individual bond energies of the Morse form, 
by counting the total number of M-M, N-N and M-N bonds in each (2 x 2) unit cell. Accordingly, for x — 
0.25, 0.5 and 0.75, we obtain the elastic contribution to the formation energy as: 

AH;j% = GVMNias) - miMias) - W^Nias) (3) 
AH^% = 8Vmn (as) - Wmm [as) - Wnn {as) (4) 
AH^% = WMN{as) - iVMM{as) - WNN{.as) (5) 

Note that Equations ([3]) and ([5]) are identical, i.e., within our model, the elastic interactions lead to a Ai?^'° 
that is symmetric about x = 0.5. It is also important to note that for bulk alloys of M and N , there are no 
terms analogous to the second and third terms on the right-hand-sides of the above equations. Due to the 
presence of the substrate, these terms have to be evaluated not at fo*^^^ or (where they would lead to a 
zero contribution) but at the substrate spacing as- As a result of this, one can expect mixing rules to be quite 
different for surface alloys than for bulk alloys; we will return to this point further below. 

In order to evaluate Eqs. ©-([S]), we need values for the Morse parameters Aq, Ai and 6, which appear in Eq. [2l 
We obtain these by computing the surface stress, a^^^^ , for each single-component monolayer as a function of 
the in-plane bondlength l^y. In principle, this could be obtained by compressing or expanding a monolayer of 
M or iV on S. However, this would make the overlayer incommensurate with the substrate, leading to a surface 
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FIG. 4: (color online) Ab initio results for AH, the formation energy per unit surface area: for each composition, the 
formation energy is plotted as a function of x, the fraction of the magnetic element. The panels on the left and right 
contain results for Fe and Co alloys respectively. Note that Pb alloys are the most stable, followed by Pt, and Ag alloys 
are the least stable. 



unit cell which is too large for practical computation. Hence, we instead compress or expand the whole slab to 
perform calculations at different l^y, and then subtract out the contribution from the substrate layers to the 

the following equation: 



total stress, so as to get the surface stress at each Ixy^ In order to carry out this procedure, we make use of 



2 



„surf (1 \ 1 („V.slab\T n\ V .bulk ^ . V .hulk^ xy 

<^ yi-xy) = 2 ^^^^ >^''- ^ ^""^ ^ ^>^^=^ 2 -"^ 3c 



(6) 



Here, cr*"''-^ is the surface stress at an intraplanar bond length l^y, for a slab with Ua atomic layers and dimension 
(which includes the vacuum) along z, and cr^a''"''' is the aa component of the "volume stress" for the slab 
(i.e., it has dimensions of force per unit area, as opposed to the surface stress, which has units of force per unit 
length) . Similarly, tr^^"'*^ is the aa-component of the volume stress for a bulk Ru cell that has been stretched 
or compressed to the same Ixy as the slab. Note that the geometrical factors in the last two terms are specific 
for the hep structure. 

The Morse parameters for an i-i bond can be extracted from the plot of a'^'^'^f versus Ixy for each single- 
component overlayer of M or N on the Ru surface. As an example, our results for the variation of surface 
stress with in-plane strain, for a monolayer of Au on Ru(OOOl), are shown in Fig.O qualitatively similar curves 
are obtained for other elements. The value of 5*^*^ or h^^ is given by the value of Ixy at which the graph 
crosses the a;-axis, while the values of and Ai are obtained by fitting the curve to an expression derived 
from a Morse potential. The values thus obtained for all seven overlayer elements are given in Table IIIII The 
value of 6" serves as a measure of the effective size of an atom i when placed on the Ru(OOOl) surface. The 
values obtained by us from calculations of surface stress are roughly consistent with estimates obtained from 
a consideration of the buckling of the surface alloys, together with a hard-sphere model. Note that for both 
the magnetic elements Af , fo*^^^ is smaller than 05, whereas for all the non-magnetic elements N considered by 
us, h^'^ is greater than as- However, the values of h are found to be different from af^i\, in some cases quite 
significantly so. This difference is due to the presence of both the surface (i.e, no neighbors above) and the 
substrate (different neighbors below). It is interesting to note that for Fe/Ru(0001), h > aj;°;|, whereas for all 
the other elements, b < al^\l This is presumably because Fe in the bulk form has the body centered cubic (bcc) 
structure with a coordination number of 8, whereas all the other elements have either the fee or hep structure 
with 12-fold coordination. As a result, only for Fe are the overlayer atoms more effectively coordinated when 
placed on a Ru surface. However, apart from such general observations, we were unable to discern any simple 
relationship connecting the values of b and a^'^ll ■ 

It remains to obtain the Morse parameters for M-N bonds. In analogy with the Lorentz-Berthelot mixing rules, 
the M-N bond parameters are assumed to have the form, 6*^^ = (6*^*^ + 6^^)/2; A^^^ = ^^mm^jva/ ^^^^ 
j^MN ^ ^^MAf ^wAr^ Table [IV] we have tabulated the values of 5*''^. These should compared to the bulk 
Ru NN spacing (= 2.74 A). It is also instructive to compare the values in Table ITVl with those in Table [H one 
finds that there is no dramatic change upon accounting for altered surface sizes. 




FIG. 5: a"'"^-' versus l^y for Au/Ru(0001): The value of a""''' at each l^y is computed by compressing or stretching the 
complete slab, and then subtracting the contribution of substrate layers. The data points are fitted by an expression 
derived from a Morse potential, to get the values of A^'", A^" and 6^". The dots represent the calculated data points, 
while the line is the Morse fit. 
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TABLE III: Values for the Morse parameters for M-M and N-N interactions, as deduced from surface stress calculations. 
The last column contains our values for the calculated nearest neighbor spacing for M or N in their bulk structures. 
These may be compared with our values for b, which is the preferred intra-atomic spacing for a monolayer of M or A*' 
on Ru(OOOl). 
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TABLE IV: Values of b''"^ , obtained by taking the average of preferred nearest neighbor spacing on the surface; these 
values should be compared to our calculated value of 2.74 Afor the NN distance on the Ru substrate. 

Our results for the elastic contribution to the formation energy, evaluated using Equations ([3]) to ([5]), are 
displayed in Fig. [51 We find that for all ten combinations considered by us, elastic interactions always favor 
mixing of the two overlayer elements, in accordance with the predictions by TersofFi^ 

At first sight, our most surprising result appears to be our finding that for both magnetic elements, the Pb 
alloys are the most stable, though upon examining Table U or Table [TVl one might think that this is unlikely. 
However, this is because for surface alloys, unlike bulk alloys, the phase-segregated forms can cost a high elastic 
energy, because of the presence of the substrate. Since pseudomorphic Pb/Ru(0001) costs a great deal in elastic 
energy, the mixed form is correspondingly favored. In order to make this argument clearer, in Fig. [71 we have 
separated out the individual contributions to the right- hand-side of Equation [H The first (M-N) term is always 
positive, while the second (M-M) and third {N-N) terms are always negative. In order for AH to be negative, 
the first term should be small (the simple mixing rule applies only to this term), while the second and third 
terms should be large in magnitude. The first term is found to follow the expectations from an elementary 
consideration of sizes (either at the bulk or at the surface): Ag and Au alloys are the most favored, followed by 
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FIG. 6: (color online) Results for AH"''^, the elastic contribution to the formation energy, as a function of fractional 
composition x, for (a) Fe and (b) Co alloys. Upon considering elastic interactions alone, Pb alloys appear to be the most 
stable against phase segregation, while Cd and Pt alloys are the least stable. 
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FIG. 7: The various contributions (from Eq. |4)) to the elastic contribution to the alloy formation energy, for x = 0.5. 
Contributions from M-N, M-M and N-N bonds are displayed separately. All the histograms have been plotted on the 
same scale, to make comparison easier. Note the very high negative contribution from Pb-Pb bonds for both Fe-Pb 
(upper panel) and Co-Pb (lower panel) alloys. 



Pt and Cd, and then Pb. It is interesting to note that both Cd and Pt alloys have roughly the same contribution 
from this first term; this is because Cd undergoes a relatively large contraction in size at the surface, relative 
to the bulk. A Co monolayer on Ru(OOOl) is relatively unhappy (i.e., the contribution to the elastic part of 
the formation energy is significant and negative), and a Pb monolayer on Ru(OOOl) is extremely unfavorable 
energetically. As a result of these two facts, elastic interactions favor the formation of Co-A^ alloys over Fe-A^ 
alloys, and lead to the high stability against phase segregation of M-Pb alloys. However, one should be cautious 
in interpreting these results, since we have made the assumption that the alloys as well as phase-segregated 
monolayers remain pseudomorphic. For the alloys, this is probably a valid assumption, since the elastic energy 
(corresponding to the first bars in Fig. [7]) is small, i.e., the stress is unlikely to be high enough to drive the 
overlayer to relax. Despite the significant elastic energy contained in a Co/Ru(0001) monolayer, it does not 
reconstruct.^"* However, for Pb/Ru(0001), the very high elastic energy makes it seem possible that this system 
might reconstruct, presumably via a network of misfit dislocations; we are not aware of any experimental 
information on this system. Thus, the high stability we obtain for M-Pb alloys may be misleading; the stability 
would be lowered if the phase segregated form were to reconstruct (since the third term in the clastic energy 
would then be decreased in magnitude). 

Finally, in Fig [51 we display our results for the chemical contribution to the formation energy, obtained by 




FIG. 8: (color online) The chemical contribution to the alloy formation energy, as a function of fractional composition 
X, for (a) Fe and (b) Co surface alloys on Ru(OOOl). These graphs have been obtained by subtracting out the elastic 
contribution from the ab initio results. 



subtracting out the data in Fig. [6] from that in Fig. [4l When only chemical interactions are considered, Fe 
alloys are more favorable than Co alloys. From this plot, we note that the stability of Ft alloys is largely due to 
the favored chemical bonds between Fc-Pt and Co-Pt ; this is consistent with the fact that these intermetallic 
systems also form bulk alloys. The Ag alloys are not stable because Co-Ag and Fe-Ag bonds cost very high 
chemical energy, which cannot be offset by elastic energy; this too is consistent with previous results . ^^'^"^i^^ 
Also note that Fe-Au bonds favor mixing, whereas Co-Au bonds cost energy, which explains the particular order 
of stability observed in the ab initio results. 



IV. SUMMARY AND CONCLUSIONS 

In this paper, we have attempted to gain an understanding into the factors governing the energetics of strain- 
stabilized surface alloys, by performing ab initio calculations on ten combinations involving a magnetic and a 
non-magnetic metal co-deposited on a Ru(OOOl) substrate. In many cases, we find the surface alloy to be stable, 
even though the constituent elements are immiscible in the bulk. 

We find that the stability (against phase-segregation) does not correlate with expectations based upon the 
simple argument that the mean atomic size should be as close as possible to the substrate lattice spacing. One 
reason for this is that though elastic interactions are an important mechanism governing stability, chemical 
interactions can also play a crucial role. In some cases, the latter are large enough to disfavor atomic-level 
mixing, even if it helps in lowering the elastic energy. A second complicating factor is that unlike for bulk 
alloys, for such strain-stabilized surface alloys, the phase segregated forms can also cost elastic energy. Thus, 
there are three factors that determine whether or not mixing takes place at the atomic level: (i) the elastic 
energy of the alloy phase, (ii) the elastic energies of the phase-segregated monolayers on the substrate, and 
(iii) chemical interactions. Because of this complicated situation, a simple criterion, analogous to the first 
Hume-Rothery rule for bulk alloys, does not seem possible for such systems. 

We have also found that effective atomic sizes on the Ru substrate are not equal to the bulk size; in some cases 
this difference is small, while in other cases it is large. Several alloys involving a magnetic and a non-magnetic 
element on a Ru(OOOl) surface are found to be stable against phase segregation; this is primarily because the 
effective size of the magnetic elements is smaller than the nearest-neighbor distance in the substrate, while that 
of the non-magnetic elements is larger, even after accounting for altered sizes at the surface. Of the systems 
we have considered, we feel that Fe-Au, Fe-Cd and Co-Cd are particularly promising candidates that would be 
worth experimental investigation. In these systems, both chemical and elastic interactions promote alloying. 
We have also found that surface alloys involving Pb and either Fe or Co appear to be very resistant to phase 
segregation; however, this conclusion is dependent on our assumption that a monolayer of Pb on Ru(OOOl) does 
not reconstruct, which may or may not be valid. 
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